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I. INTRODICTION 

Conventional nethods usea to represent underwater sound 
transmission have depicted transmission loss (a3) as a func- 
mon Of range in either tabular or graphical form. Since 
the interest was in determining quantities such as range of 
bottom ktounce, convergence zones, and VepEooabt lity 9.0L 
detection, these methods proved very satisfactory. However 
when it became necessary tc compare various underwater 
acoustic propagation ioss models witn each other or actual 
in-situ experiments, they provided the analyst with less 
than adeéeguate insight into the causes for the okserved 
errors in the various nodels outputs. Furthermore, interest 
has been generated in analyzing the receivei signal for 
information from the sound source which cannot easily be 
induced from transmission loss as a function oi range. 

One method, the Wavenumber Technigue (WT), described by 
Meme DiNapoli of NUSC [ke£f. 1] and applied ty Rickard Lauer 
of NORDA [Eef. 2], provides more information for model 
comparisons and nay have tne capability of Getermining sound 


source depth and range from the receiver in certain cases. 


DiNapoii used an analysis of acoustic propagation in the 
wavenunber domain ag€an il1nternediate step in getting the 
transmission loss curve in the rast Field Progran BEE). 


Lauer studied the wavenumber domain information trom the FFP 
Meemconcluded that source localization might be vossible. 
Peo. otamey ; Ref. 3] investigated the potential for deter- 
Mining tke source depth by using the Brock Split-Step Fast 
Bem@rter Transrorm (SSFFT) [Ref. 4] parabolic 2guation model 
With tne "Wavenumber Technigue" (WI) model. His preliminary 
investigation involved the use of isospeed and ASTREX sound 
Beeed prozil2?s with a fuliy absorbing botton, varied 


source/receiver combinations, and aultipie frecuencies. 
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The present investigation uses two acoustic farabolic 
equation models for comparison with theoretical Llovd nirror 
depictions, and the effects of observed inconsistencies in 
predicted transmission loss curves on the calcuiated results 
from application of the WT model. An analysis of the sensi- 
tivity in determining source depth based on model output is 
attempted to ascertain possible use in naval operations. 

Momracrlitate a Straight-forward analysis of the results 
and evaluation of the WT technigue for the cases studied, 
only a selected aumber of the variouS computer runs aré 
macluded. Jnly those curves which best illustrate the 


effects revealed py this research aré poresentéd. 
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II. WAVENUMBER TECHNIQUE (WT) 


A. BACKGROUND (LLOYD MIRROR) 


The WI can be elucidated with the help of the classical 
Lloyd mirror effect [Ref. 5] which describes the interaction 
between direct path and surface reflected sound signals [ron 
the same continuous wave (CW) source. Figure 2.1 illus- 


trates the geometry involved in the Lloyd mirror effect. 


©. 
ZA 
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na 2 a | 
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Figure 2.1 Lloyd Mirror Geometry. 
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The ccmplex pressure (P) can be expressed by 


1 -jkr -} s 
a KR ai 


P2A| 4 


eo 


when ©, and r, are related by 


a S 
fi 
(egn 2.1) 


rE = Pace r= Seen: Keane) 


with the assumption r,;~ Lr, the pressure amplitude (p) can be 


approximated by 


Ry = Sind ieee (egn 2.3) 


The interaction between the direct path and surface 
reflected waves produces constructive and destructive inter- 
ference which is manifested as peaks and nulls for the in 
phase and out of phase conditions respectively. aCe: zie. 2 
waS produced by using eguation 2.3. tae nul Pe locations can 


be oktained by the relationship. 


2A (Ans 
-_—-— —— {(egn 2.4) 
Dp R sin R 4 


when equation 2.3 1s expressed in the horizontal wavenumber 
domain (k-,), the details of the resultant functional depen- 


dency on horizontal wavenumber yield useful information. 
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B. PHYSICAL PROCESS 


The WT is 
{with range 
Spectral density as a 
AT 
a guadrature 
signal the 
demodulation process can best be 
6 j. 


dependent amplitude 
mito the 
wavenumber. The use or the 
ment would reqjuire 


to attain complex 


meme rgure 2.3 { Ref. 


2cos (ut) 







-2sin (Gat) 


Figure 2.3 


When acoustic models 
the 


eliminated 


reguirement to perform a 


since the 


i= 
Oo? 


a process by which the 


ee 


filter 


LP 
filter 


are used to simuiate 


complex 


complex pressure wave 
and phase) is transformed 
EUnet non Ob “the horizontal 
in the operational environ- 
demodulation of the source 
The 


defined py the illustration 


pressure. quadrature 


Acos|(w - We) t + x] 





A sin|(w ~GA)t+ x| 





Quadrature Demodulation of a Signal. 


the environment, 


quadrature demodulation is 


pressure 1s directly 





accessible. Pr ctUd hey meeGOustlC Pdbabolic equation 
models provide the complex envelope function (U) which must 
be multiplied by the Hankel function in order to obtain the 


complex pressure. The acoustic wave 1S written in the forn 
P=UHS’(k r)eJut (eqn 2.5) 


where the reference wavenumker is defined by 


e 





kKo=-e;5 (eqn 2.6) 
and 

C=O ee) (eqn 2.7) 
fmemtne Solution to the appropriate parabolic equation. The 


complex pressure is corrected for volune attenuation and 
Fourier transformed to attain tne spectral density. A oor 
Gisplaying the results graphically with the horizontal wave- 
humber on the x-axis and the normalized spectral density on 
the yraxis is constructed and the spacings between the nulls 
measured. Figure 2.4 {from Ref. 2] 15 an example of the NT 
for a source and receiver in the same type of water mass. 

Since the Lloyd mirror field in tne wavenumber domain is 
given by [Ref. 2] 


chee oe ae (eqn 2.8) 


ny 





an alternative formulation [Ref. 2], a#hich produces evenly 
Seaced nulls, is obtained by converting the horizontal wave- 


number to beta 


B=\V ki -k? am 


By using beta instead of the aorizontal wavenuaber, the 
distance, celta beta, CDetween any two nulls can be used to 


ascertain the source depth by 


sp =f (eqn 2.10) 


aS illustrated in figure 2. 5. 

The application for implying tne nzethod of images, in 
isospeed cases, to describe the erfects of a perfectly fiat 
pressure release boundary can be justified by inspection. 
Therefore the parabolic equation can be used to produce the 
Meoyda's mirror effect and surface decoupling 1S not an 
issue, cf. Chapter 4 of Brekhovskikh [ Ref. 7]. 

When the WT is implemented, care must be exercised to 
ensure that the complex pressure contains an adequate number 
of points to describe the shortest periodicity in the pres- 
Sure field. Gils ewiles beclude the possibility of aliasing 
in the wavenuaber donain. The real and imaginary elements 
ef the pressure signal must be modified so that the signal 
begins and ends with zero values. ims modetication is 
necessary because tne pressure signal must represent a 
repetitive oscillaticn for the transform. The complex pres- 
Sure array is zero-filled beyond tne data in order to 
generate an array which has the length of a power of 2 for 
the Fast Fourier Transform (FFT). The norizontal wavenumber 


increment is generated within tne computer code by 
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III. ACOUSTIC MODE 


ANALYSIS 


a a aE 


Ae PARABOLIC EQUATICN 


The elliptical wave eguation can be approximated bya 
parabolic equation (PE) when it is assumed that the envelope 
function varies slowly with range. A detailed nathematical 
description of the PE from an acoustic point of view is 
presented by DeSanto [Ref. 8 and 9] and a simpler but less 
general description can be found in Coppens { Ref. 10]. The 


parabolic equation has the fora 





2 
é Ye + 2k (ZY) + (k?-K2)U =O (ean 3-1) 
Z 


where it is assumed that the pressure has the form of equa- 


mean 3.2, 


Dee geez Se) (eqn 3.2) 


S(t) represents the primary radial dependence of the field 
in terms of an outward propagating cylindrical wave of the 
form [Rerf. 4] 


aA 
Spe Kc) (cane ces) 


If it 1s assumed that the range of interest is many wave- 
lengths from the source; then the asymptotic form of the 


Hankel function (equation 3.4) can be used [Ref. 11]. 


ee 





(1)_,/_2 9 j(ker-/4) 2 (eqn 3.4) 
Hs =V>e> e? , Kor => 1 ean 


Two PE models, the Brock version of the Split-Sterp Fast 
Fourier Transform (SSFFT) and the Jeager version of the 
Implicit Finite Difference (IFD) were investigated for 
compatibility with the WT in determination of source depth. 


Bae SSFFT 


1. ransmission Loss Comparison with Lloyd Mirro 





To analyze sound waves in the wavenumber domain, we 
needed a model to provide the pressure as a function of 
range. Initially the SSFFT model [Ref. 4] was selected. 
This model generates successive values for U as a function 


of range with the helr of the algoritha 
U(r+ear,z)=e jarke(n*-1)/2 


(eqn 3.5) 


FFT | jark? 
e JATK72Ko FET U(r 7z)) 


This model was a natural choice since it was the model used 
Mae ne preliminary study dealing with the determination of 
the source depth by the WT [Ref. 3]. However, 1t soon 
became apparent that this model's inherent weaknesses would 
require the consideration of another nodel if ocean botton 
interactions were to be studied. These weaknesses, and 
their impact upon the WT, will be discussed shortly. 

The SSFFT 1S a range dependent acoustic wave model 
which, for this analysis, will be operated ina range inde- 
pendent manner. In other words, Only one sound speed 
profile will be used, the water mass will be assumed homoge- 
neous, and the bottom flat. TOmracliiztate the study of the 
SSFFT model, the source code of the nodel was modified to 


generate the variables required for follow-on programs to 


ae 





Beenmsrorm the pressure information and tnen display the 
results graphically. ihecomryeOrE tac SOUrCe Code listing ior 
each program is provided in appendices A, 3B, and C. These 


programs were then linked by Job Control Language ({JCL 


UV) 
O 


that the WT graphics were an automatically executed. 

A test cun of tne SSFFIT was made for 100 az, 500 on 
source and receiver depths, 0.019 km range step, fuel 
Maser ring botton, 2000 an water depth, and an isospeed 
profile at 1590 n/sec. These values were chosen so that a 
Bomeparison colld be made with the thaoretical Lloyd mirror 
transmission ioss (figure 2.2). Tne output was displayed as 
meeransnisSion loss curve (figure 3.1). PEO Go ulation 2.4, 
mmerleccation 9f null number 1 (R1) snould be 33.33 xa and 
null numrer 2 (RK2) should be at 16.67 xn. The SSrFI placed 
feeeac a2proximately 38.3 km and R2 at 17.1 kn. Sance th 


(D 


model kept the frequency, source depta, and receiver Cepth 
constant, then this equates to a rererence sound speei of 
1305 and 1462 m/sec respectively. hunts dood s ene, that. a 

the range decreases the error is reducei which indicates the 
Possibility of a problem in the range step section o£ the 


model which leads tec the introduction of a systematic and 


hth 


Me@ereasing deviation of the calculated sound field from that 


Meedicted by the classicai Lloyd's mirror interterence 


n 

pattern. ta2 “washing out" of the interrerence pattern for 
Bapeges less tran about 2 km isa result of the particular 
MeeeoxiMation or a point source used to initiate the 
program, which is tc Le expected. Ratner regular ztluctua-~ 
tlons develop bevond about 30 xa; BROUmGRe GCONCtEEY OF-tne 
case, it 1s plausible that this arises from interference 
Witn bottom reflected signals. 
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Another test was made uSing the same inputs as above excendt 
that the source and receiver were aoved to a depth of 1500 2 
feaGure 3.2). Since the half beamwiith is 20 degrees and 
the bottom is fully absorbing, tne Snortest range at which a 
Mjerace ceflection ccuid occur 1S approximately 8.2 ka; 
beyond this range, surrace reflections should be spaced in 
accordaace with equation 2.4. fae nulls beyord &.2 kn do 
occur at the correct ranges; however, there are siow mnodula- 
tions of the overall transmission loss signal. These Siow 
Mmea@ibations are another indication that the bottom is tot 
Tuliy aksorbing as had been specified as input to the 


progran. 


ZG 
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2- 4&2 Analysis 


Piomiimecitout from, the SSFrT model is in Eig 
units, waich iS a probiem when comparing different nodeis. 
However, for this investigation the problem waS minimal. As 
expressed by equation 2.6, the reference wavenumber for 100 
hz with a reference sounad speed of 4921.26 ft/sec (1500 
Mesec) is 0.1277 1/ft. After correction of tne reference 
wavenumber for a range increment of 62.3 it and the 2048 
points used to represent the wavenunter spectrun, the 
Maximum beta value should be Q.1243 1/ft which agrees 
reasonably well with the maximum peta viiues in figures 3.3, 
or, and 3.5. SOlebary TtONtneOn, che oOULlS acre not equasly 
Spaced and the expected delta beta value of V9.00192 1/it 
(equation 2.1)) is net represented by any of the null spac- 
gS . Stamey [Ref. 3], observed tne same Non-unitora 
spacing or reta nulls. The reasoa why the source deoth is 
not represented by any or the nuli Spacinys may be tracei to 
tne improper null placement in the raage donaii. Tae Lack 
of equal spacing between the nulls must be attribut¢a to 
something else, possibly the pressure wave was igproperly 
bepresented as a continuous Signal or numerical errors exist 
Meenin the source ccde. Whatever the reason for this 
difference, it must be resolved bezore the SSFFI can be used 


With tne KT to determine the Source depth. 
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IFD 


-_ 


the IFD nodel (ref. 12 ) whicn was gastalles and 


sasted if September 4983 at NPS { Ref. 13] was the next nodeL 
nosen as 4 sogund source for the Wi. ie decision Was made 
ecause mye tFrD "a> capable or properly handling the oceah 
orton interactions and was the only otaet Noctst ne == nodet 
n resilience at NPS- Neo fy of the source code listing Ss 
provided in appendix ee 

The IFD also presented prooiens which wile ne 
discussed further. The same initial test ran ysea wita the 
SSEFT Was ysed with the ILFD- Since tne TED cequz-c= gore 
potton information, addit ional testing was oerfiormed in 
order to obtain a fap, absorbing botton. EaguLe pee LS 
trates tne Ponder. o> seen by the Trp model. 

The densities in the watet masSe sediment, ine ae 
ficial attenuat 10b ay coe ee set to 3 constant Bae end 
the attenuation in the sediment was adjasted antil the el1a- 
gation of botcton interrere nce was onpserved. Figure sa 2s 
an exaaple oz tne IFD +ransmissioa loss output for a 300 2m 
source/receiver depth, 4 water Geptn of 2009 
surface of the artificial attenuation layer at 2000 G, a @ 
lower pressure release surface at 1900 n- te botton atten 


yation used to produce figure 3.7 was 9.00 16 gnywaveleng: be 





SO Ce 


water 


sediment 


“~upper surface of 
artificial attenuation layer 


nce inser oat 


pressure release surfaces 





Figure 3.6 IFD Environment. 
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faemocation of the nulls in ‘figure 3.7 coincide with those 


Wa 


predicted by equaticn 2.4 if tne higa frequency noise i 
paeperea out by either a running average process or low pass 
Gmgatal filter. A careful examination of the transmission 
Weeeecurve (figure 3.8) reveals the real problen, which is 
the long wavelength noise near tae wavelength of interest. 
“re Long wavelength noise will cause erroneous spikes on the 
Meepeots after being Fourier transformed. Goisimted anbal= 
Meee revealed that the a9ise could be eliminated if the 
thicknesses of the artificial attenuation iayer and thre 
bottom pressure rcelease surface are chosen so that they are 
at least twice the water depth. Piece lOles OL a. saat l 
vertical grid step will reduce the thickness reguired for 
the pressure release surface, but this version of the IFD 
Beogram is linited to 5000 vertical grid points due to soft- 
ware restrictions. MEGuecman outs a thansilssion Loss curve 
obtained from a 500 m source/receiver depth, a 2000 m water 
mptn, a400) mn upper levei artificial attenuation layer, 
and an 3000 m bottom pressure release surface. Comparison 


MemeereuresS 3.1) and 3.9 wito the classical Lioyd's mirror 


Medals @ission loss (figure 2.z) snows that tne Shape of the 
Subve beyond 30 km is greatly improved with the IFD. Tae 
Number or vertical grid steps was 5990. ahen the IFD is 


Meed in a shallcw water situation and at 25 hz, tne noise 
Mees NOt aprear t> be present. Figure 3.10 is a plot of the 
Beaiomission loss for a Signal at 25 hz, Witn the source/ 
Beceiver deptn at 50 mn, the bottom sloping up from 350 to 50 
M, and an artificial attenuation layer and pottom pressure 
Memedse Surface located at 750 and 1000 a, respectively. 
Pe noCise whith was so evident in Ligure 3.7 is absent fron 
the Shallow s#«ater case. Producing tre neise free curves 
requires excessive computer tine except for shallow water 


Cases. 
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py 


2. HT Analysis 


To reduce the amount of the pressure information and 
yet preserve the detaiis, the IFD tests were conducted 
between 10 and 40 km in range. The reference wavenunber for 
Mmreecest input with the IFD was 0.4189 1/fm which, after 
Sgerecting ror a range increment of 7.5 m with 4096 points, 
produced a beta maximum of 0.4188 1/m. Pgures. S241, 3.12, 
and 3.12 indicate agreement with the beta maximun ard a 
@eea beta of 0.00623 1/4 can be found among the null spac- 
ings. When che vertical grid step 1s reduced in size, the 
mmeput signal contains 2ore information of finer detail. 
Mie finer detailed information will, if noise is present, 
Meodguce more clutter on the WT depiction (figure 3.13). en 
Numerous WI depictions a peak was observed at the Dinianun 
and maximum wavenumbers. Staney believed that the right and 
left intensity maxima corresponded to bean elevation angies 
of 0 and 30 degrees, respectively. He rturther stated that 


Beewlett fiaxina also may be reiated to algorithmic inaccura- 


emes {Ref. 3]. The "U-shaped phenomenon" observed during 
Stamey*s investigaticn or the SSFFI, is present with tne IFD 
moe the nulis are not equally spaced. 
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ftmee ceason for the “U-shaped pnenomenon" appears to Oe 





celated to noise in the envelope function field Which is 
i amplified when tke Hankel function is inserted to 


Ompain the conplex pressure. Using a low pass filter ana 


lincreasing the vertical grid step size appear to reduce the 


/problen. iment aetic DULIS OGCeir at 33.30 and 33.39 


‘km when a single null should exist at 33.33 km. This condi- 


i 


|tion oor noise produces the most devastating erfects by 





| masking the desired density spectrua. A fartow band pass 


filter would probably solve this problem in the short tern 


but eventually the IFD will require further study to elini- 


nate the noise. When the shallow case is used in the WT 
(figure 3.14), the "U-Shaped phenomenon" is adDsent and the 
null locations are more distinct. However, i1n the shallow 


water case, the pottom interaction produces high frequency 
interference but the null spacings occur hear tne expected 


value for the source depth (0.0623 1/m). 
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HORIZONTAL WAVENUMBER 


1 meter Vertical Grid Step. 


Shallow Water WT, 


IFD 


Figure 3.14 





™WO parabolic equation models, tne SSFFT and IFD, were 
Meed tO predict sound fieids for comparison with the Licya's 


mirror interference fattern in the rang? domain and then in 


tJ 


the wavenumber domain. Mir ecie Eange €demain, the SSFr 
improperly handled the location orf pressure nulls ana 
Seesgolayed bottom interrerence. While these weaknesses in 
the SSFrTI were absent in the IrD, the inadvertent and dele- 
terious insertion of noise was extensive. In tae wavenunber 
Memain, the noise in the IFD was prohibitive in ail but the 
déep sediment and shallow water cases. tio: Deby prioduced 
results which were ccmparable with the SSFFI wnen the noise 
WaS not a factor. The range step and tae botton attenuation 


megorithms of the SSFFT should be reviewed in order to 


correct those problems descrired acove. Pieaene FLED, the 
memes propiem will require turther investigation. te toe 
Peoplem in the IiFD 1s not corrected, ir SemOge taal be 
severely limited. With the guidance provided by Lauer 


(Ref. 2], the wavenumber domain can provide tne detailed 
information needed Soueg Geek. correct inconsistencies 
opserved petween acoustic models. THE Tez Ore, the FT should 
be investigated further to ensure taat any sortware weax- 
hesses are eliminated and a thorough understazainsg o£ the 
process is assured. 

ExperimentS in the operational efvironment Should obe 
meammed to test the NT with actual sound sources. The shis- 
Foard experiment could be easily acconplisned since the WT 
Can be 1aipiemented on micro-computers if the vressure infor- 
Mrmemon 1S provided from an outside source. fhe wT, if it is 
gOing to be considered as a metnod to determine source 


depth, appears to have three major snortconmings: 
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1. A clean CW Signal with pnase information is essen 
tliat, wWhteh will regquite quadrature demodulation and 
extensive filtering. 

2. Prediction podels are reguired which will provide 
equaliy spaced beta nulls so tnat software can be 
developed to automate the process. 

3. High resolution processing eguipment is necessary 
because the change 1n beta spacing for shallow sound 


sources is small (figure 2.5). 
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